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Abstract 

In c = 1 units the product (mass x radius) for the neutron and 
the proton is about 4.7/1 assuming their radii equal to lfm. We show 
that the corresponding products for the Dirac neutral and charged 
membrane coincide and are equal 1.6h. 

1 Introduction 

In 1962 Dirac, trying to explain the value of the muon mass, con- 
sidered an idea that the electron could be modeled by a conducting, 
elastic membrane of spherical topology [1] (at the same time he also 
considered a neutral membrane in the presence of the gravitational 
field [2] - the hnite size of the electron was originally considered much 
earlier by Lorentz, Abraham, Bucherer and Langevin, see [3] and ref- 
erences therein) . Assuming the Lagrangian L = Lem + Lmem. , where 
Lem is a usual term for the electromagnetic held and Lmem. is given 
by the membrane world-volume, it can be shown that for the spheri- 
cally symmetric case with the radial coordinate p, the Hamiltonian of 
the model (in c = 1 units) is given by 

H, = J-h^dj + uj^p^ + ^^ (1) 



2p 

where oj/in is the membrane tension and e is the electric charge. Us- 
ing the Bohr-Sommerfeld quantization method, one hnds that the hrst 
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excitation of the membrane corresponds to the energy ~ 53me , where 
nie is the mass of the electron. Dirac notes that in order to get a 
closer value to the experimental ~ 207me one would presumably 
have to introduce spin into the theory by pertorming the square root 
appearing in Hr- While this might be true, in this paper we focus 
on another possible application of Dirac's membrane i.e. the effective 
description of hadrons. We find a surprisingly good agreement be- 
tween the product (mass x radius) for charged/netral membrane and 
the experimental values for the proton/neutron respectively. 



Among all closed (compact, without the boundary) objects coupled 
to the electromagnetic held in D = 4 spacetime dimensions, mem- 
branes play a special role ensuring that A^^ is nowhere singular. The 
reason for this lying in the fact that such object split into two 
disconnected regions: the interior and the exterior. Equations of mo- 
tion for membranes are nonlinear hence quite complicated and unlike 
in the case of points or strings very little is known about the exact 
solutions (for some examples see [1|). As for the quantum theory, the 
existing method of quantizing the bosonic membrane [5] is highly non- 
trivial resulting in a certain quantum mechanical model with matrix 
degrees of freedom. When the size of a matrix is taken to inhnity the 
precise description of the quantum membrane is obtained. In spite of 
these complications recently there has been a considerable progress in 
understanding the bosonic membrane theory [6]. 

Since closed membranes in provide a natural split of space into 
the interior and the exterior there exists a proffered curvilinear system 

in spacetime and a function /(x) such that the equation /(x) = 
describes a membrane and equations /(x) > 0, /(x) < describe a 
region outside or inside the membrane, respectively. It is convenient 
to fix the curvilinear system such that x^ = f{x) and choose = 
x^ =: T, = = x^ for the three variables o"", a = 0, 1, 2 - the 

internal parametrization of the membrane word-volume. The action 
for a membrane coupled to the electromagnetic held considered by 
Dirac [V\ is (in c = 1 units) 



2 Dirac's model 



S = Sem + Si 



'mem. i 



(2) 




X 
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Smem = / Mdx^dx^dx^ 

4vr 7^.1=0 

where g^'^ is the metric corresponding to the curvihnear system x^ 
(concretely Dirac takes the induced metric g^y = d^y^dyy\, A = 
0,1,2,3, where y^ are rectihnear and orthogonal), J = ^ — det g^^i, 
and M = J\/—g^^- The couphng is due to the factor J appearing 
in both Sem and Smem.- Yarying ([2]) with respect to y^ one arrives 
at the equations of motion which for the sphericahy symmetric case 
x^ = r — p, x^ = 0, x^ = (so that the surtace is given by x^ = 0) 
give 

d P ^ 2 

dt yr3-2 2a;p4 pyr— 2- 

The balance between the repulsive electromagnetic forces and the at- 
tractive ones, due to the positive membrane tension, is when p = 
hence = /4lo where a is the radius of the electron. On the other 
hand the total energy of a system at rest E = e^/2/> + /3p^ (minimal in 
the equihbrium provided /3 = w) is equal to both 3e^/4a and mg when 
p = a. Theretore one concludes that a = 3e^/4me = 0.75re where 
Te is the classical electron radius, ~ 2.8fm (the value a = 2.1fm 
is of course not realistic as it is bigger then the charge radius of the 
proton 0.87fm [7] - according to Dehmelt the radius of the electron 
could be of order lO^^^lm [8j). Considering small oscillation about the 
equilibrium one hnds the corresponding frequency to be \/6/a hence 
the energy of one quantum would be hu = \/6h/a = 448me. 

To improve the analysis one proceeds to the hamiltonian formal- 
ism. While the details of this are quite involved the hnal answer for the 
spherically symmetric case turns out to be particularly simple given 
by ^ ■ The complications are due to the choice of the coordinates s'^ 
hence the loss of the explicit diffeomorphism invariance of the action. 
There exists a generally covariant formulation [9] in which one con- 
siders, in addition to Sem and Smem. written for an arbitrary ^'^ and 
(7°, an extra term proportional to e^da^^^A^ where is the current- 
charge density on the membrane. One then shows that all Dirac's 
hndings, in particular the hamiltonian ([1]), can be obtained in a less 
elaborate way keeping the explicit diffeomorphism invariance. 

Using the Bohr-Sommerfeld quantization condition one obtains an 
approximation which can be written as 

mpirac ^ 1 / 32V^r(7/4) \ i 
me 3 V ar(l/4) ^ 
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where a = /hc is the fine structure constant, obtaining 



52.4, 83.1, 109.0, 132.0, ... 

me 

for consecutive values of n. 

This approximation can be improved as observed in [TO] by noting 
that the in the a — )• hmit the coulomb term in ([T|) acts like an inlinite 
wall and hence one should apply different boundary conditions in the 
Bohr-Sommerfeld quantization procedure. As a result one obtains 

mimpToved ^ ^g .^^ ^Q2.9^ 126.1, ... 

me 

which agree with the numerical values also found in [10]. 



2.1 Further investigation 

Let us rewrite the hamiltonian Hr in ([T|) in dimensionless variable 

X = p/a 

where we used u) = -^, a = 

4a'' ' 4me 

While hnding an analytic expression for the square root of the 
positive differential operator H is in general dithcult the calculation 
of the matrix elements of the square root y/H is possible using the 
standard procedure by diagonalizing H and taking ^/H = U^^^/ĔU 
where U is s.t. H = U^^EU, E = diag{En)- Applying numerical 
techniques (for details see the Appendix - we use a different technique 
compared to |10] ) we hnd that this prescription gives 

^^^w 43.6, 76.3, 103.0, 126.5, ... 

nie 

which are consistent with the results of [lOj . 

In order to find the radius of the hrst two excitations (let us denote 
the corresponding wavefunctions by -i/^^ for Ei = 43.6me and ipr for 
E2 = 76.3me) we use their probability densities (Figure 1). It lollows 
that \'>pfi{p)\'^ has a maximum for p ~ 6.7a while [^/'^(p)!^ has two 
maxima for p ^ 3.8a and p ^ 11. 9a so that the probability density 
\^pr{p)\'^ will be viewed in as two concentric membranes. 

It is somewhat interesting to consider the generalization of Dirac 
action to the case of two (or more) concentric membranes. We will 
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Figure 1: Probability densities and l-^T-^a;)^. 

assume that the total charge e is on the exterior membrane. If the 
tension of the inner, neutral shell is k(jj, k > the corresponding 
hamiltonian of this system would be 

Hrp = V-d? + k^oj^r^ + + a;V^ + ^ (5) 

The eigenvalues of this hamiltonian are hfted compared to Hr- One 
can choose k so that its hrst excitation is equal to the mass of the 
muon (this takes place for k ~ 63) however the second excitation is 
not even close to the mass of the taon. Of course one could consider 
three concentric membranes and fit the additional tension so that the 
mass of the taon appears but such model would be just a fit. 



3 Reinterpretation 

Perhaps the result = 43.6me could be improved by introducing 
spin in some clever way into the theory. However it seems that this 
membrane model works surprisingly well if we identity charged mem- 
brane with the proton and the neutral membrane with the neutron. 

Let us start with the case of the neutral membrane. The discrete- 
ness of the hamiltonian ([1]) is not due to the Coulomb term. Clas- 
sically, without the electromagnetic field the membrane will collapse 
(i.e. p = after a finite time - p{t) can be obtained from the equa- 
tions of motion which imply pp = 2(/9^ — 1) solved by the Jacobi sine 
tunction) but the operator 

Hneut. = ^-h^dj + n^p^ = {h^^f'''^-dl+y^ (6) 

where y = p/{h/il.)^'^ is dimensionless, itselt is discrete (since it is a 
square root of a discrete, positive dehnite operator) hence quantum 
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mechanically a free membrane will develop bound sates. The natural 
unit of energy is now {Jt?^})^'^ where the tension is not specitied. 
The hamiltonian Hneut. has parity even and odd eigenstates. The 
parity even states have the maximum of the probabihty density for 
/9 = which corresponds to a membrane with radius - a point. 
However the parity odd states have the maximum for p > for which 
the eigenvalues are 



{h?Qy/'^ 



The probability density of the first excitation mo ~ 1.95{h^Q)^/^ has 
a maximum at po ~ 0.82{h/Q)^/^ (Figure 2) hence mopo ~ 1.6/1. 



I<^m„(y)l 




Figure 2: First excitation for a neutral membrane. 

It is quite amusing that for the case of the neutron we have m„p„ ~ 
4:.7h (in c = 1 units, we took pn = lfm) in rough agreement with the 
result for the neutral membrane. Considering the fact that we do not 
mention quarks, gluons and the spin, this naive picture of a neutron 
as a neutral membrane works surprisingly well. 

We can now do the analogous analysis for the charged membrane. 
The hamiltonian is the same as for the case of the electron only written 
in different units (cp. (jl])) 



Hcharged = {Ah^il/a)'/' -dl + + ^ j 

(where, z = p/{ah/An)^/^ , n can be related to the classical radius at 
the equilibrium) . Due to the Coulomb term we must assume il^{0) = 
for the wavefunctions hence all the excitations have nonzero radius. 
The eigenvalues are 

mcharged _ ^.239, 0.418, 0.564, 0.693, ... 



(4^2Q/c,)l/3 
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and the lowest energy state mo ~ 0.239(47i^O/a)^/^ is peaked at po w 
6.75(0/1/4^1)^/^ (same as for ip^^, see Figure 1) implying mopo ~ 1.6h. 



4 Conclusions 



It is not a surprise that effective description of hadrons in terms of bag- 
like models will essentially give results consistent with experimental 
data. However what is striking about the Dirac membrane model is 
its huge simphcity compared to the actual processes taking place - 
no quarks, gluons, spin - one wonders why this model works at ah? 
Moreover, the value of the only dimensional parameter of the theory, 
the tension Q, was nowhere used as it cancels out in the calculation. It 
seems that, at least for the purposes of calculating the product (mass 
X radius), the spin can be neglected while the strong forces replaced 
by a 2 + 1 dimensional held theory - signaling a sort of holographic 
[12j behavior. 

Acknowledgment I thank M. Kuzniak and P. O. Mazur for discus- 
sions and the correspondence as well as the Swedish Research Council 
and KTH for support. 



we use the cutoff method [11] which consists of calculating the matrix 
elements of in some orthonormal basis (on [0,cxd) in our case), 
truncate the inhnite matrix and then numerically diagonalize it. The 
spectra and the eigenvectors of the truncated matrices converge very 
quickly to their exact counterparts. 

An important step in this method is the choice of the convenient 
basis (a common choice is the basis in the Fock space - not suitable 
here). Due to the l/x part in the regularity of the wavefunction 
il){x) at the origin implies that il){x) ~ x' for / > 1. For this reason 
we choose the orthonormal basis en{x), n = 1, 2, 3, . . . as 



Appendix 



In order to obtain numerically the eigenvalues of 




en{x) 



^n{n + 1) 



1 



xL\'_^{x)e "^, / en{x)em{x)dx = 
Jo 




nm 
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where Ln^ {x) are generalized Laguerre polynomials. In this basis the 
explicit representation of operators K = —^i, W = x^ and V = Ij^ 
can be obtained with 

At this point we introduce a cutofF Njnax, consider a Snite matrix 
h^nm"'''^ = + Wnm, n,m < Nmax ^'iid then numerically find its 
eigenvalues E^^""""^ and the matrix J/^^""*^) s.t. 

j^{Nmax) _ Jj(Nmax) T -^(Nma^) Jj{Nmax) ^(Nmax) _ (Hdg^E^^"^'^^^ ) 

The square root of (^g^j^ obtained and the matrix rep- 

resentation for the overall operator is 

Jj{Nmax) ^ u{Nmax) T ^i^rnax) jj(^N^^^) _^ y(Wmax)_ 



Thc spcctrum of H: 



quickly convcrgcs to thc exact values (Figurc 



3). To obtain the wavefunctions corresponding to ^;^^™"^) the 
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Figure 3: The conyergence of the eigenyalues of H-, 



(Nmax) 
X 



eigenvectors 



of iJ^^-»-) i.e. 



N„ 



k=l 



The convergence of '(p^"""'\x) to the exact i/jni^) is governed by the 
corresponding convergence of the eigenvectors. In practice the cutoff= 
10 already gives very accurate approximation for first excitations. 
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To find the eigenvalues of the parity even states of ([6]) we use the 
same numerical method but with help of a different, more convenient 
basis 

fn{x) = =H2n~i{x)e~'' /^, / fn{x)fm{x)dx = 6nm 

where Hn{x) are Hermite polynomials. 
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